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Abstract 

We study the single transverse-spin asymmetry for inclusive light-hadron productions in 
the proton-proton collision, p^p — )■ hX {h = it, K, 77), for the RHIC kinematics based on the 
twist-3 mechanism in the collinear factorization. The analysis includes all the contributions 
from the soft-gluon pole and the soft-fermion pole for the twist-3 quark-gluon correlation 
functions in the transversely polarized proton. After discussing the flavor decomposition 
and the P^-dependence of the asymmetry obtained in the previous analysis at the center- 
of-mass energy = 62.4 and 200 GeV, we will give a prediction for the asymmetry at 
y/S = 500 GeV and also for the 77-meson production. We found slightly smaller asymmetry 
at = 500 GeV for 7r^'° and compared with those at the lower energies. The 
asymmetry for the 77- meson turned out to be significantly larger than that for vr''. 



1 introduction 



In the last decades, large single transverse-spin asymmetries (SSA) in (semi-)inclusive reac- 
tions have been receiving much attention in high-energy spin physics (for a recent review, 
see [1]). Since the large SSA can not occur in the conventional framework for the high- 
energy process based on the parton model and perturbative QCD, understanding the origin 
of SSA provides us with a new opportunity to reveal the nucleon structure and the QCD 
dynamics. One of the most conspicuous SSA is the one measured in the inclusive single- 
hadron production in the proton-proton collision: 

pHp,s^)+p(p')^h(Pk) + x, (1) 

where S± is the transverse-spin vector of the polarized proton and h = n, K and f] etc. The 
SSA for this process is characterized hy = {a'^ — a^) / {a^ + a^) = Aa/a, where a'^{a^) is 
the cross section obtained with the proton polarized along the spin vector S± {—S±). The 
FNAL-E704 collaboration reported the first data for An for p^p — ?> ttX and p^p -> nX at 
the center-of-mass energy ^ = 20 GeV, which is as large as 30 % in the forward direction 
of the polarized proton or antiproton [2, 3, 4, 5]. RHIC at Brookhaven National Laboratory 
also reported a similar magnitude of An at VS =62.4, 200 GeV [6, 7, 8, 9, 10, 11]. 

In our recent paper [12] (hereafter referred to as KKIO), we have presented a numerical 
analysis of the RHIC An data for p^p hX {h = vr, K) at y/S = 62.4 and 200 GeV in 
the framework of the collincar factorization. Since the unpolarized cross section for this 
process has been well-described by the next-to-leading-order (NLO) QCD analysis in the 
coUinear factorization [13], the analysis oi An in the framework can be taken as a guide to 
clarify the origin of the asymmetry. In the coUinear factorization, SSA appears as a twist-3 
observable and is described in terms of multiparton correlation functions [14, 15, 16]. (For 
the detailed formalism of the twist-3 calculation demonstrating the gauge invariance and 
the factorization property, see [16, 17, 18].) In KKIO, we focused on the contribution 
from the quark-gluon correlation functions Gp{xi,X2) and Gp{xi,X2) with quark-flavor 
a {a — u, d, s, u, d, s) in the transversely polarized nucleon. They are defined from the 
light-cone correlation functions of the form ~ ['^"■gF-^^ip"') in the nucleon, and xi and 
X2 — Xi represent, respectively, the light-cone momentum fractions carried by the quark 
and gluon lines coming out of the nucleon. These functions satisfy the symmetry property 
G'^p{xi,X2) = Gp{x2,Xi) and Gp{xi,X2) = —Gp{x2,Xi). (For the definition and the basic 
properties of these twist-3 functions, see KKIO and [16, 19].) In the twist-3 mechanism for 
SSA, the contribution from these two functions to the single-spin-dcpcndent cross section 
Act occurs from a pole part of an internal propagator in the hard part, which reflects the 
naively T-odd nature of SSA. Such poles are classifled into the soft-gluon-pole (SGP) and 
the soft-fermion-pole (SFP), which flx the momentum fractions at Xi — X2 and = (i —1 
or 2), respectively. Because of the symmetry property, Gp does not contribute through SGP. 
Therefore the cross section can be written as [20, 21, 22, 23, 24] 

Aa = E (g^x, x) - x^^^^) ® fix') ® D%z) a^Xc 

a,b,c \ / 

+ Y: (G^(0, x) + G^p(0, x)) ® fix') ® D%z) ® (2) 

a,b,c 



where f^{x') is the twist-2 unpolarized parton distribution in the unpolarized nucleon 
and D'^{z) is the twist-2 unpolarized fragmentation function for the final hadron with 
the parton labels h and c (6, c = q,q,g)- The scale dependence of each nonperturbative 
function is implicit. As shown in (2), the SGP contribution appears in the combina- 
tion of Gp{x,x) — a;— [22, 23], while the SFP contribution appears in the form 
G%{0,x) + G^(0,x) [24]. In KKIO, we have performed a fit of the RHIC data for the 
pion and kaon productions based on Eq. (2) and found that the existing data have been well 
reproduced by the combination of the two effects. Although the main contribution is from 
the SGP contribution, the SFP contribution also plays an important role which can not be 
substituted by the SGP one. Since we could not get a good fit by the SGP contribution 
only, we shall focus on FIT 1 of KKIO which includes both SGP and SFP contributions 
below. This analysis used the GRV98 for the unpolarized parton distribution [25] and the 
DSS fragmentation function [26] for tt and K. 

The purpose of writing this paper is twofold. The first is to present more details of the 
analysis of FIT 1 in KKIO [12]. In particular, we explain the role of the SGP and SFP 
functions for each quark fiavor leading to the obtained pattern of Ajv for the pions and 
kaons. We also show the dependence of on the final hadron's transverse momentum 
Pt in wider Pt region and explain the origin of the characteristic P^- dependence. The 
second purpose is to give a prediction for A^ for pions and kaons at = 500 GcV and 
also for the rj meson at y/S = 200 GeV, which are being measured at RHIC. Confrontation 
of these predictions with the forthcoming RHIC data is critically important to test our 
parametrization and to clarify the mechanism of the observed A^^s. 

During the present study, we found that the cross section formula for the SGP contri- 
bution used in KKIO (Eq.(6) of [12]) has an error in the overall sign. (We have been using 
the convention for the covariant derivative = — igA^ and the e-tensor €0123 = 1, 
but Eq.(6) of [12] is not consistent with this convention. ) Since the SGP and the SFP 
functions are the independent functions, the correction of this sign error can be done by 
simply reversing the sign of the extracted SGP functions without affecting any results for 
the obtained asymmetries. For completeness we will show in Sec. 2 the corrected SGP func- 
tion for each quark and antiquark flavor together with the SFP functions. In the meantime, 
the sign error in [22] for the extracted SGP function was also reported in a recent paper 
[27]. By our present correction, our SGP cross section formula becomes consistent with 
what is claimed in [27]. 

The remainder of this paper is organized as follows: In Sec. 2, we present the flavor 
decomposition of the obtained A^ for the pions and kaons, and discuss the role of each 
function in the asymmetry. In Sec. 3, we show the PT-dcpcndcncc oi A^ for 7r° in the 
wider Pt range and clarify the origin of its characteristic behavior consistent with the 
RHIC data [8]. In Sec. 4, we will give a prediction of Ajv's at \/S — 500 GeV using 
the parametrization of FIT 1 of KKIO, anticipating its measurement at RHIC in the near 
future. In Sec. 5, we present a prediction for A^^ for the rj meson at y/S = 200 GeV, which 
reproduces the observed tendency of the preliminary RHIC data [10]. Section 6 is devoted 
to a brief summary. 
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Figure 1: Decomposition of at Xi;' > obtained in FIT 1 of KKIO [12] into the SGP 
and SFP contributions (left panels) which are decomposed further into each quark-flavor 
component (right panels). The solid and dotted lines show the contributions from the SGP 
and SFP components, respectively. The dashed line in the left panels shows the total A^^. 
Data points from RHIC [8, 9] are also shown. 



2 Flavor structure of the asymmetry 

We begin our discussion by the flavor structure of the asymmetry for 7r^'° and K^. In Fig. 1, 
we show the decomposition oi Ajq into the SGP and SFP contributions for each quark and 
antiquark flavor. The left panels show the decomposition of into the total SGP and SFP 
contributions, while the right panels show the contribution from each quark and antiquark 
flavor both for the SGP and SFP contributions. To understand the pattern of the flavor 
structure, we show the SGP and SFP functions obtained in FIT 1 of KKIO for each quark 
and antiquark flavor in Fig. 2, by correcting the sign error mentioned at the end of the 
introduction. First, as seen in the top two panels of Fig. 1, the major contributions for tt^ 
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Figure 2: The SGP function and the SFP function G^(0,x) + G^(0,x) at the 

scale = 2.4 GeV^ for each quark flavor obtained in FIT 1 of KKIO in comparison to the 
unpolarized quark distribution fa{x) [scaled by factor 1/10]. 



come from the SGP components for the favored quark flavors: Ajy for vr"*" is dominated by 
the SGP contribution for the u and d quarks, and Aj^j for vr" is dominated by that for the 
d and u quarks. Note here the significance of the u and d components, which results from 
the large SGP functions for these "sea" flavors in the nucleon at large x as shown in Fig. 2. 
The net SFP effect turned out to be small for vr^ in the region of xp > 0.4 for which the 
RHIC data are available. If one looks into more details of the SFP contribution, the SFP 
effect brings some net effect for tt~ , while it is much smaller for vr^. This feature of the 
SFP contribution can also be seen for K"^. (See discussions below.) The above feature of 
the SGP dominance can also be seen in Ajsf for vr*^: All the SGP contributions from u, d, u 
and d bring large contributions to A^^ in the large xp region. The net SFP effect can be 
found in the small xp region for ix^: It brings some effect at xp < 0.4 which cancels the 
SGP contribution, making the An small. This feature can be understood if one remembers 
that the partonic hard cross section for SFP is much larger than that for SGP in many 
relevant channels, although the derivative of the SGP function enhances its contribution. 

The twist-3 functions for the sea quark flavors in the nucleon become more important 
for the kaons. In particular, the observed large A^ for K~ should come from the sea quarks 
due to its flavor structure, K~ ~ ms. As seen in Fig. 1, the largest contribution to A^ for 
and K~ comes from the SGP functions for the s and s quarks, respectively, and for 
even the SGP contribution from the w-quark is significantly smaller than that from the 
s-quark. This is due to the much larger strangeness component for the kaon in the DSS 
fragmentation function as well as the sizable magnitude of the SGP functions for the s and 
s quarks in the nucleon. Combined with the result for the pions, one may say that A^- for 
the pions and kaons reflects the flavor structure of the fragmentation functions rather than 
that of the SGP functions. This is in contrast to the unpolarized cross section where the 
contribution from the valence quark flavors in the nucleon dominate in the forward rapidity 
region. 
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Figure 3: Decomposition of the w-quark SGP and SFP contribution to Aat at \fS = 62 A 
GeV shown in Fig. 1 into the quark and gluon-fragmentation channel. 



The contribution of the SFP effect shows a very different pattern from the SGP contri- 
bution. As shown in Fig. 1, the net SFP effect survives for K~ and 7r~ at xp > 0.4, while 
it is negligible for and 7r+. This net effect for the former is caused by the SFP func- 
tion for the u-quark flavor through the gluon-fragmentation channel in combination with 
the large gluon component in the DSS fragmentation function, in particular, for the kaon. 
(The SFP function for the d-quark is also large and partly cancels the u-quark contribu- 
tion. However, the w-quark function has a longer range than the ci-quark function, and thus 
the former effect is larger.) As shown in [24], the SFP partonic hard cross section in the 
gluon-fragmentation channel is extremely large and has the opposite sign compared with 
that in the quark fragmentation channel. For and 7r+, the w-quark SFP contribution in 
the gluon-fragmentation channel is mostly canceled by that in the w-quark-fragmentation 
channel for which K'^ and vr^ have a large fragmentation function as a favored flavor, while 
the former contribution is not canceled by the latter for K~ and vr". To see this feature, 
we showed in Fig. 3 the decomposition of the SGP and SFP contribution for the u-quark 
into the quark-fragmentation channel and the gluon-fragmentation channel. As seen from 
this figure, the u-quark SFP contribution via the gluon-fragmentation channel survive for 
K~ and vr", while they are canceled by that via the quark-fragmentation channel for 
and 7r+. 

Summarizing the flavor structure of A^, the SGP contribution for the favored quark 
flavors is the dominant one, while the SFP contribution for the valence u and d quarks is 
also significant via the gluon-fragmentation channel. This SFP effect survive for K~ and 
7r~ but is canceled by that in the quark- fragmentation channel for and 7r+. 



3 Py-dependence 

The P^-dependence of provides an important test for the mechanism of the asymmetry. 
In Ref. [12], we have found that our predictions for the 7r° production shows a good agree- 
ment with the RHIC-STAR data at = 200 GeV. In Fig. 4, we show the P^-dependence 
of Aj^ for 7r° in the wider range of Pr at xp = 0, 0.28, 0.37 and 0.50 together with the 
RHIC-PHENIX data at xp = [7]. As seen in the figure, calculated A^ at xp = stays 
zero in the whole Pp region, which is in agreement with the PHENIX data. 

A conspicuous feature of the P^-dependence is that the calculated A^- once increases up 
to Pp ~ a few GeV and then decreases slowly, which is in agreement with the STAR data. 
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Figure 4: P^- dependence of A]^ at = 200 GeV and four fixed xp together with the 
PHENIX data at = [7]. 



This behavior is in conflict with our naive expectation that should behave like 1/Pt 
as a function of Pt by reflecting its twist-3 nature. We have found that the contribution 
from the gluon-fragmentation channels plays an important role to cause this peculiar Pt- 
dependence of A^- It brings a negative contribution to the positive polarized cross section 
Act which is the numerator of Ajsi, while it brings a positive one to the unpolarized cross 
section a. With the de Florian-Sassot-Stratmann (DSS) fragmentation function used which 
has a large gluon component, their contributions decrease very rapidly as the scale /i = Pt 
becomes higher. Correspondingly, Aa decreases relatively slowly around Py ~ 1 — 3 GeV 
compared with the denominator, even though the twist-3 cross section accompanies the 
factor 1/Pt- Such behavior is absent in the previous study by Kouvaris et al [22] which 
used Kretzer's fragmentation function [28] whose gluon fragmentation function is small. 
For comparison we calculated A^ by fixing the scale of the distribution and fragmentation 
functions at yU = 1 GeV (instead of = Pt) and found that Aj^ decreases monotonously as 
Pt increases and diverges toward P^ — J- 0. 

Another interesting point in Fig. 4 is that An does not decrease as fast as 1/Pr in 
the high Pt region. This is because of the coexistence of the two effects proportional 
to M^Pt / {—T) and M^Pt / {—U) in the twist-3 asymmetry where T = {p — P^)"^ and 
U = {p' — PhY are the Mandelstam variables: While the former effect is much larger at 
small Pti the latter also brings significant contribution at large Pt- 

4 Aat at = 500 GeV 

Since the measurement of A^r is being extended to = 500 GeV at RHIC, we present a 
prediction for An aX this energy, using the SGP and SFP functions obtained from the data 
at = 200 and 62.4 GeV. 

In Fig. 5, we show the result oi An for pions and kaons as a function of for the 
pseudorapidity t] = 3.6 and rj = 4.0. For the pions, the general trend of the asymmetry is 
the same as those at lower energies, while the magnitude of the asymmetry becomes smaller 
at each xp. This is because a larger value of y/S corresponds to larger Pp for a fixed xp 
and thus the polarized cross section is more power suppressed at higher \/S. 
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Figure 5: A]\f for pions and kaons at = 500 GeV and the pseudorapidity rj = 3.6 (left) 
and 7] = 4.0 (right). 

A]sf for kaons turned out to become larger than those for charged pions in the large 
xp region, in particular, for the K~ meson. This larger for is mainly due to the 
faster decrease of the unpolarized cross section for the production compared with the 
vr^ production with increasing y/S from 62.4 to 500 GeV. This decrease of the unpolarized 
cross section becomes even more manifest for the K~ production, since the unpolarized cross 
section in the forward region is dominated by the contribution from the valence flavors in 
the nucleon and the sea contribution is more suppressed at higher scale fi = Pt a.t large 
xp- Therefore A^- for K~ becomes much larger than that for K~^. We found that A^- for 
K~ overshoots one in the large xp region. In our present analysis, the constraint on the 
SGP and SFP functions for the s and s quarks comes only from a very limited number of 
the kaon data with the assumption that they are the same for the s and s quarks. With a 
larger number of data points, this unfavorable feature should be corrected. Nevertheless, 
we expect the relation A^ > A^'^ > A]^ at y/S = 500 GeV for the reason stated above. 

Finally we show in Fig. 6 the Py-dependence oi A^ for vr" at = 500 GeV. As in 
the case of 75 = 200 GeV, An once increases toward Pt — a few GeV and then decreases 
slowly. 

5 An for r/ meson 

The first data of An for the 77- meson production was reported by the FNAL-E704 collab- 
oration for the pp and pp collisions [5]. They found a clear indication that for the pp 
collision is positive and substantially larger than A]^ . Recently the RHIC-STAR collab- 
oration reported a preliminary data for A]^ at = 200 GeV, which also shows A^ is 
larger than A]^ in the large xp region by about factor two or more [10]. Here we present an 
estimate for based on Eq. (2) using the SGP and SFP functions obtained in KKIO. In 
this scenario, the difference between A]^ and A]^ comes from the difference between the dif- 
ference in the fragmentation functions. As the fragmentation function for the //-meson, we 
adopt the recent Aidala- Ellinghaus- Sassot- Seele- Stratmann parametrization [29] which 
is obtained by NLO QCD analysis of the existing data, assuming the equal contributions 
from all quark and antiquark flavors as = = = = = Dg. Although we 
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Figure 6: Pr-dependence of A]^ for VS = 500 GeV at xp =0, 0.28, 0.37 and 0.5. 

used the leading-order (LO) formula for the twist-3 cross section (2), we use the above only 

available NLO parametrization for the //-fragmentation function in this first study of A^. 

1 

Figure 7 shows the calculated A]^ as a function of Xi? in comparison with at ^/S = 
200 GeV and rj = 3.65. In the large xp region, A]^ is approximately 1.5 times larger than 
A]^ . This is consistent with the observed tendency of A]^ at RHIC, although its magnitude 
still seems insufficient to explain the RHIC-STAR data [10]. As in the case of Ajf , the 
SGP contribution dominates A^. To see the origin of the difference between Ajf and AJ^ , 
we have shown in Fig. 8 the decomposition of the SGP contribution to A]^ and A]^ into 
the contribution from each fragmentation channel. From this figure, one sees that the 
difference between the two mesons comes from the large contribution in the strangeness 
fragmentation channel for the //-meson. This large strangeness contribution is caused by 
the combination of the large strangeness component in the ?7-meson fragmentation function 
and the large SGP functions for s and s which were responsible for A^ . Therefore in 
the present analysis, there is a strong correlation between the large Aj^ for the kaon and 
the ?7-meson. As another possible origin of the difference between and Aj^ , the twist-3 
fragmentation function may contribute [30], for which the complete cross section formula 
is not available yet. This effect may explain the remaining difference between the two 
asymmetries. For the complete clarification of the asymmetries, one needs to perform a 
global analysis of a greater variety of data using the complete cross section formula, which 
is beyond the scope of this study. 

6 Summary 

In this paper, we have studied the SSA for the inclusive pion, kaon and r^-meson productions 
in the pp collision at typical kinematical regime of RHIC, using the SGP and SFP functions 
for the twist-3 quark-gluon correlation functions obtained in our previous analysis. We first 
discussed the role of the SGP and SFP functions for each quark and anti-quark fiavor. It 



^ Since the LO and NLO parametrizations for tt in [26] are not very different from each other, we expect 
the use of the NLO fragmentation function for rj would not cause much error. 





Figure 8: Decomposition of A]^ and A]!l into each fragmentation channel at ^/S = 200 GeV 
and the pseudorapidity t] = 3.65. 

turned out that the SGP contributions from both "valence" and sea flavors in the nucleon 
play an important role to reproduce the observed pattern of the asymmetries, and that 
the main contribution to An for each meson reflects more the flavor content in the DSS 
fragmentation function. Our analysis reproduces the observed Py-dependence of RHIC An 
data for 7r°. We found that the large contribution from the gluon- fragmentation channel 
both in the unpolarized and polarized cross sections plays a crucial role to cause the peculiar 
behavior of A^ around Pt — a. few GeV. 

We also presented a prediction oi A^ for the kaons and pions at ^ = 500 GeV, which 
is being measured at RHIC. We found a slightly smaller asymmetry for the pion compared 
with those at y/S = 62.4 and 200 GeV. For the kaons, Ajy was found to become larger than 
those for the pions in the large XF-iegion at = 500 GeV, especially A^ for K~ becomes 



much larger because of the small unpolarized cross section. An for the 77-meson turned out 
to be larger than that for 7r° by about factor 1.5, which partially explains the tendency of the 
RHIC data. This larger originates from the large strangeness fragmentation function for 
the 7] meson combined with the sizable strangeness SGP function in the polarized nucleon. 

Present analysis is based on the assumption that the whole A^ for p^p hX {h = tt, K 
and r]) originates from the twist-3 quark-gluon correlation functions in the polarized nucleon. 
To determine these functions, the SSA data on the Drell-Yan process and the direct-photon 
production in the pp-collision is extremely useful, since there is no fragmentation ambiguity 
associated with the final hadron in these processes. For the complete understanding of 
SSAs, a global analysis of variety of processes based on the complete formula is needed. 
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